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stabilized by shell effect against the instability due to repulsive component of nuclear force, 

Tij" ' analogous to superheavy elements similarly stabilized against Coulomb instability. Here we 

Sp^ ' discuss the stability and fission decay properties of such nuclei in the ^^^Th region and show 

^^ ' that they are stable against a and fission decay and have /3-decay life time of several tens of 

seconds. In particular, the ^^^Th nucleus has a low fission barrier and unusally large barrier 

width. This makes it an ideal thermally fissile nucleus, if formed by means of a thermal neutron, 

O I like other known nuclei such as ^^^U, ^^^U, ^^^Pu in this actinide region. It shows a new mode 

r^ I of fast fission decay, which may be termed as multifragmentation fission, in which in addition to 

two heavy fragments large number of scission neutrons are simultaneously produced. Its likely 

L^ ■ synthesis during the r-process nucleosynthesis will have important bearing on steller evolution, 

^ • and here in the laboratory, it has great potential in energy production. 
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The superheavy elements in Z = 114 — 126 region Ql 12 El which are vulnerable 
to spontaneous fission [U 13 E] due to the infiuence of repulsive Coulomb force 
of protons, have been predicted to be stable due to shell effect. Due to this effect, the 
periodic table has extended upto Z=116 [T by now with promise to extend further along 
the valley of stabilitiy, thus elongating the stability peninsula El E]- Another effect 
that was shown by us earlier is the shell stabilisation of highly neutron rich superheavy 
isotopes oi Z = 62,78 and 90 nuclei ^2]- It is known that the nucleon-nucleon force 
itself has a repulsive component (triplet-triplet, singlet-singlet) whose contribution pro- 
gressively increases with the increase of neutron number making the nucleus unstable. 
The shell effect can stabilize this instability giving rise to new magic numbers in the 
vicinity of neutron drip-line. This is a complementary parallel effect which should widen 
the stability peninsula. In fact, in an extensive study involving three different methods, 
infinite nuclear matter (INM) model ^3 El El El? relativistic mean field (RMF) the- 
ory 11711111012011211 and Strutinsky shell correction calculation ^21, it was shown [T2I 
that islands of stability around new magic numbers N=164, Z=90; N=150, Z=78; and 
N=100, Z=62 will exist giving rise to superheavy isotopes ^^"^Th, ^^'^Pt and ^^^Sm. The 
N/Z ratios of these nuclei are 1.82, 1.92 and 1.62 respectively, in contrast to the value 
of 1.54 for the doubly closed shell superheavy element "^^^Xu^ which is quite similar to 
that of ^''^Pb and ^^^U in the valley of stability. This result is depicted here in Figure 1. 
It may be recalled that a new frontier in 1980 's has opened up with the production of ra- 
dioactive ion beams with the prospect of synthesis of about 5000 nuclei in the laboratory, 
and thereby, extending the nuclear peninsula upto drip-line region. The ultra neutron 
rich nuclei are expected to show totally new properties and will probably be of immense 
utility to mankind. Here we investigate the properties of a representative ultra-neutron- 
rich nucleus ^^^Th which is a superheavy isotope of thorium. It is known that ^^^Th is 
a fertile material whose fission decay properties are well known. We show that ^^"^Th is 
stable against a— decay and fission decay and has /?~decay half-life of tens of seconds. 
Most interestingly, ^^^Th has unusual fission properties with low fission barrier of 3.57 
MeV, but a very large barrier-width, which makes it infinitily stable against spontaneous 



fission but fissionable if formed by means of a tliermal neutron. It will undergo a new 
mode of fission decay in wliicli, in addition to two fragments, large number of surplus 
neutrons will be instantaneously produced. This may be termed as multi-fragmentation 
fission. Due to highly stable character, this ^^^Th nucleus will have important impli- 
cations for r— process nuclear synthesis in stellar evolution and also great potential for 
energy production in the laboratory. 

Alpha and Beta decay of Th isotopes 

The usual magic nuclei with nucleon numbers 2, 8, 20, 50, 82 and 126 have their 
special properties, in particular extra stability, due to the close shell, and they occur 
along the central line in the valley of stability. However, as one goes across the valley 
towards the neutron drip line, as is shown in Figure 1, new shell closures can arise 
giving rise to new magic nuclei, and new islands of stability around them. The search 
and discovery of such nuclei is on the frontier of nuclear physics research today. It will 
be interesting to study the decay properties of the ultra neutron rich heavy nuclei. It is 
well known that, in general, heavy nuclei with A > 150 are unstable against a— decay 
while for lighter nuclei this process is improbable. Infact the superheavy nuclei ^^^Xii4, 
^^^Xii4 show such decay, which seems to be the case for all heavy nuclei along the line 
of stability. The heavy neutron rich nuclei are expected to show deviation from this 
feature. In Figure 2a, the a— decay half-lives of different Th isotopes calculated using 
the predictions of the masses in INM |24j and Howard-Moller (HM)[25| mass models 
and RMF theory [21] are shown. The calculation of a— decay half-life T^ has been 
done using the widely used standard formula logT^/ sec. = (aZ + b)Q~^^'^ + (cZ + d), 
where a = 1.66175, b = -8.5166, c = -0.20228 and d = -33.9069, due to Viola and 
Seaborg ^7\. The value of Qa is estimated by using the standard relation Qa{Z,N) = 
BE{Z -2,N-2) + BE{2,2)-BE{Z,N), where BE is the binding energy. In general, all 
the isotopes are quite stable against a— decay with decay hfe-time exceeding 10^°years 
for all the three mass models. The RMF and INM show a prominent peak around 



N = 160 ~ 168 as mark of extra stability for isotopes around N=164. Beyond that 
the Qa values become negative prohibiting such decay altogether. Thus ^^^Th can be 
considered to be stable against a— decay. 

The neutron-rich nuclei are generally f3~ active. They will successively emit /^"particles 
until they reach the line of /3— stability in the valley, a parallel scenerio to the superheavy 
elements which attend stability by successive a— emission (see Figure 1). The /?— decay 
half-lives T^ are estimated by using the standard formula due to Seeger, Fowler and 
Clayton J2H| ^/3 = ^^^^wd '"^ sec, where A =number of states of the daughter nucleus 







within 1 MeV of ground state nucleus x exp{— A/290) with A being its mass number and 
Wo = BE{Z + l,A) - BE{Z, A) + 1.29 MeV. The half-lives of Th isotopes so obtained 
are presented in Figure 2b, which stop falling rapidly from N=164 with some flat peaks 
thereafter indicative of extra stability, and the values are in the realm of tens of seconds. 
Thus ^^"^Th can be considered to be a fairly stable nucleus for laboratory studies. 

Fission decay of ^^^Th 

The study of the fission decay properties of ^^^Th is of special importance because 
the lighter stable isotope ^^^Th is known to be a fertile nucleus. The question how the 
ultra-neutron-rich superheavy isotopes for the actinide nuclei will decay by fission has 
not been addressed before. It is here one may expect some exotic new phenomena. 
The fission barrier determines the fission decay properties of the nucleus. The nucleus 
can overcome the barrier and undergo decay when excited by an external agency like 
an energetic neutron incident upon it. The neutron separation energy Sn provides a 
measure of excitation for the fission decay process. We have plotted in Figure 3 the 
neutron separation energy Sn and fission barrier Bj of the even Th isotopes as a function 
of neutron number. The black squares and the green stars are the S'„ values of the HM 
and INM models. They are connected by lines to guide the eye. The closed red circles 
represent the values of fission barrier Bf taken from Howard and MoUer [23] . It is well 
known that if Sn < Bj, then the nucleus can not undergo thermal neutron fission. The 



fission tliresliold Enm = Bf — Sn has to be overcome by impiging with an energetic 
neutron, and thereby the nucleus will undergo fission. Hence, if Sn > Bf, then thermal 
neutron (with practically zero energy) can cause fission. The thermal fission becomes 
a very attractive possibility for energy generation. It is interesting to see in Figure 3 
that, although normal Th (^'^^Th) is not thermally fissile, its isotopes around N=164 
in the range N=162— 170 are thermally fissile, which has important implication and 
consequences. We have as yet only a few generally known nuclei which are fissionable 
with thermal neutron, out of which only one i.e., ^ss-jj jg naturally occuring with an 
abundance of 0.71% in natural uranium, and the others like ^ss-y ^^^^^ 239p^ ^^^ ^^ 
artificially produced. Of course they have long half-lives of 1.6 x 10^ and 2.1 x 10^ 
years respectively, but ^^^Th as will be shown aposteriori is permanently stable against 
spontaneous fission. 

Fission barrier profile of ^^^Th 

We now consider the nature of the fission decay mode of ^^"^Th, which is primarily 
governed by the profile of the fission barrier. The height and width of the fission barrier 
which is supposed to be parabolic in nature have to be obtained. We have followed an 
empirical method to get the width of the barrier from the systematics of the known 
experimental fission half-lives, and extrapolate them to extremely neutron rich region 
of interest to us. The fission half-lives can be calculated in a simplistic manner as 
Ti/2 = ln2/np, where n is the number of barrier assault by the decaying fragment, 
related to the width hu by nh = huj/2'K. And p is the penetrability of the barrier given 
by p = [1 + exp{2Ti B f / huoY^ . Taking the values of Bf from reference [221, and using the 
experimental ri/2 [29j, we get the values of huj, the systematics of which so obtained are 
shown in the plot of huo versus Bf in Figure 4 for various actinide nuclei. It is indeed very 
revealing as well as interesting that the plot shows a linear behaviour with progressively 
increasing slope with the increase of proton number of the elements from 90 to 96. For 
the next element Cf with Z=98, the linear behavior gets fuzzy. However, the mean 



follows the trend with a higher inclination. The trend is more conspicuously restored for 
Fm with Z=100. This deviation correlates well with the fission mass yield systematics, 
where considerable deviation from the standard well-defined two peaks occur for Fm 
isotopes. It may be mentioned that such systematics of Bf versus hu has been analysed 
for the first time here. The width huj for any isotope with calculated fission barrier may 
be obtained by extrapolation of the linear graph. For ^^^Th with a fission barrier of 3.57 
MeV, we obtained the value of hu close to zero. Now we can construct the parabolic 
barrier of base width Ar somewhat schematically, using the value of hu following the 
relation l/hu = (PV/dr"^, where V is the potential energy. The barrier so obtained is 
presented as the blue curve in Figure 5, which is an extremely flat and wide barrier. 
For comparison, we have presented in the same figure the fission barrier of ^^^Th as 
red curve constructed with hu = 0.415 MeV obtained from the experimental half-life 
|29j . and Bf = 7.40 MeV from Howard and Moller [23j. The unusually large fission 
width derived above in case of ^^''Th is quite understandable and in confirmity with 
expectation. It is not necessarily true that the fission barrier falls with the increase of 
neutron number in an isotopic chain, it may even rise jHO] • So the wide barrier obtained 
here is rather specific to ^^^Th, which may be related to its shell closure structure. 

Our experience from the fission studies of uranium, has shown that the excess neu- 
trons mainly inhabit the neck region when the nucleus deforms on its journey to fission, 
which is related to the width of the barrier. Since in the case of ^^^Th, 22 excess neutrons 
participate compared to that in ^^^Th, a larger neck must be formed giving rise to a flat 
and wider barrier. This feature is also in tune with the fact that in the limiting case of no 
Coulomb force, the barrier will disappear reminiscent of infinite width and zero height. 
This interesting feature of the barrier of ^^^Th with low height but large width points 
out a new feature where, the nucleus is infinitely stable against spontaneous fission. But 
with a slight deposition of energy with a thermal neutron, it will undergo fission. Thus 
the fission half-life of spontaneous fission is nearly infinite. 

New mode of fission decay 



We next examine the decay mode of fission of a nucleus, {A, Z) decaying to {Ai, Zi) 
and {A2,Z2), solely guided by the Q— value systematics of the process, where Q is 
defined as Q^{A,Z) = BE{Ai,Zi) + BE{A2,Z2) - BE{A,Z). As is well known, the 
major driving force for the decay is the Q— value of the reaction. The probability of 
fragment mass yield in a given channel is directly related to the Q— value. In Figure 6, 
we have plotted the Q— values of the binary decay into two fragments Ai and ^2, as a 
function of the mass number of the Ai fragment, for all the relevant elements with even 
values of Z, starting from 46 to 58. The complimentary fragment {Z2, A2) is thereby 
fixed. Since the yield falls rapidly with the decrease in Q— value for an element, we have 
only shown the distribution of Q— values lying above 90% of the highest values. For 
the sake of comparison, the Q— value distributions for both ^^^Th and ^^^Th are shown 
in Figure 6. The black and red curves pertain to ^^"^Th and blue and green to ^^^Th 
calculated using the predicted masses in INM and RMF models respectively. Since ^^^Th 
is an ultra-neutron-rich isotope, during its scission many neutrons will be left free as they 
can not be bound to either of the fragments which may be already on the drip-line. In 
the figure, vertical lines mark the drip-line of the respective elements. As expected, the 
drip-line falls nearly in the middle of the distribution in case of ^^"^Th, indicating that 
all the isotopes lying to the right of it will be unstable against instantaneous neutron 
emission from the fragments at scission. These drip-line neutrons will be simultaneously 
produced along with the newly formed fragments. In the usual fission process, neutrons 
are produced by the emission from the fragments after they are fully accelerated. But in 
the present case of ^^"^Th, a certain number of neutrons will be simultaneously produced, 
signalling a new mode of fission decay which may be termed as multi-fragmentation 
fission. An order of magnitude of these prompt multi-fragmentation neutrons can be 
estimated from the mass yield plot for ^^^Th as shown in Figure 6, which turns out to 
be about 2 to 3 neutrons per fission. These are the additional neutrons apart from the 
normal multiplicity of neutrons emitted from the fragments. In case of neutron-induced 
fission of normal ^^^Th and ^35^ nuclei, the fission neutron multiplicities are of the order 
2.3 to 2.5 |HHt I34j . This number is, therefore, nearly doubled in case of ^^^Th fission, 
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which will have important implications on the energetics of the fission process. 

Although, the beta decay life-time of ^^^Th is few tens of seconds it is much larger 
than the nucleon decay life time, which is of the order of 10"^^ seconds and therefore it 
will have a great implication in the r— process nuclear synthesis and consequently stellar 
evolution. We may also conjecture that in the laboratory, such nuclei can be of immense 
potential as energy source. 

Conclusions 

In conclusion, we have studied the decay properties of Th isotopes, in particular 
^^■^Th, a superheavy isotope of thorium which has been predicted earlier to be a doubly 
close shell nuclei with N=164, and Z=90 stabilised by shell effect against repulsive 
nuclear force instability. It is stable against a— decay and has /3— half-life of several tens 
of seconds. The fission decay properties of such neutron-rich nuclei in particular in the 
actinide region have not been addressed before, which has important bearing on the 
RIB programmes in many laboratories in the world seeking the synthesis of about 5000 
nuclei in the exploration upto the drip-lines. On the basis of the systematics of neutron 
separation energies and fission barriers, it has been shown to have a rare property of being 
thermally fissile putting it into the rank of other such generally known nuclei like ^^^U, 
^^^U and ^^^Pu. Its fission barrier profile determined empirically using the systematics 
of experimentally known so far in this actinide region fission half-lives and barrier of the 
actinide nuclei, shows it to be very fiat and wide. Such feature is in confirmity with the 
expected long neck due to the excess number of neutrons to the tune of 22 over that of 
the ^^^Th, presumably supported by its shell closure property. Coupled to this, its small 
barrier of 3.57 MeV, makes it extremely stable against spontaneous fission, but highly 
vulnerable to thermal neutron fission, a quite unique property indeed. A new mode 
of fission decay named multifragmentation fission is predicted for this nucleus, where 
in addition to two heavy fragments, 2 to 3 scission neutrons will be simultaneously 
produced. This is in addition to the normal multiplicity of neutron emitted by the 



fragments which are of the order of 2.3 to 2.5. Thus the doubhng of the neutron emission 
per fission will have strong implications as for the r— process nucleosynthesis in the steller 
evolution. It also presents an attractive possibility as a source of energy production in 
the laboratory. 
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Figure 1: The new islands of stability identified around the magic numbers N=100, 
Z=62; N=150, Z=78 and N=164, Z=90 have been shown by the blue stars in the nu- 
clear chart. The red area represents the nuclei with experimentally known binding 
energies j23]. The neutron and proton drip- lines shown as solid curves are the predic- 
tions of infinite nuclear matter model. The superheavy elements occur along the line of 
/3— stability and is shown by red circles, whereas the superheavy isotopes occur close to 
the neutron drip-line. The superheavy elements undergo a— decay and the superheavy 
isotopes /?— decay to attain higher stability. The shell effect stabilizes the superheavy 
elements against Coulomb instability while it stabilizes the superheavy isotopes against 
repulsive nuclear force instability. 
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Figure 2: (a) a— decay half-life for Th isotopes using the empirical formula logTa = 
{aZ + 6)(5^^/^ + (cZ + d) of Viola and Seaborg J2Z|. The Qq,— value calculated by using 
the mass prediction of INM, RMF and HM models. The a— decay half-lives obtained 
by all the three models show similar trend. All the three models show a peak around 
A^ = 160 — 170 before Qa becomes negative with increase of neutron number forbidding 
the process. From the y-axis it is clear that the Th isotopes are extremely stable against 
alpha— decay, (b) The /3— decay half life for Th isotopes using the formula |SI1 EHI 
Tf3 = ^^^^^ti '"^ sec. is calculated, where A =number of states of the daughter nucleus 
within 1 MeV of ground state nucleus x exp{— A/290) with A being its mass number 
and Wo = BE{Z + 1,A) - BE{Z,A) + 1.29 MeV. The value of the /3-decay half- 
lives decreases almost exponentially with increase of neutron number. ^^^Th has an 
appreciable life-time of tens of seconds. 
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Figure 3: Fission barrier Bf, and binding energy of the last neutron Sin as a function of 
mass number A for Th isotopes. The Bf are taken from Howard and MoUer [2^1 and 5'i„ 
are taken from Refs. [SHj, |21] and [S2] for HM, INM and droplet model, respectively. The 
neutron separation energy predicted by all the three models predict a general pattern. 
In all the three cases the fission barrier lies above the single particle energy, except for 
the region A^ = 162 — 170, which will be vulnerable to thermal neutron fission. 
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Figure 4: The barrier width Tiuo as function of fission barrier Bf. The Bf values and the 
experimental half-lives taken from Howard and Moller [JSj and Ref. [211; respectively 
have been used to extract the barrier widths of all the known nuclei showing fission 
decay in the actinide region. This Figure shows that for each element hu and Bf show 
a linear relationship through a straight line behavior depicted in the figure. The star on 
the line for Thorium in figure, denotes the value of huj to be 0.05 MeV for its Bf = 3.57 
MeV. This linear relationship has been used to obtain hu for unknown isotopes by 
extrapolation shown by straight line. 
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Figure 5: The schematic curve is constructed using fission barrier height Bf and width 
hu to represent the shape of the parabohc fission barrier using the expression y = —ax^. 
Here a = hu, y = cPV/dr'^ and x = Ar the base width of the parabola, where V 
is the potential energy. This figure shows that smaller the value of hu larger is the 
width. The normal nuclei like 230,232rj]^ ^^ 235,238^ have hu ^ 0.4 - O.SMeV. In the 
figure, the red curve represents the fission barrier of ^^^Th and the blue one of ^^^Th. 
For the neutron rich superheavy isotope of ^^"^Th, the width is fiattened unexpectably 
and the height decreases considerably as whown in the blue curve. This fiattening of 
fission barrier makes the nucleus stable against spontaneous fission decay, because of 
decreasing penetrability. On the other hand, due to the decrease of fission barrier, with 
a minute inducement by a thermal neutron, fission decay will occur. 




Figure 6: Q-value distribution given by Q^ = BE{Ai, Zi) + BE{A2, Z2) - BE{A, Z) 
for 232rp]^ g^j^^ ^^^Th as a function of Ai fragment in the binary decay A ^f Ai + A2. The 
binding energy used for calculation of Q— value is taken from Ref. ^| and [2^] for INM 
and RMF respectively. The fission yield decreases drastically with increase or decrease 
of mass number of given element {Ai^Zi). Therefore, we have shown the distribution 
in the range 90 to 100% of the peak value in each case. The vertical line marks the 
neutron drip-line for the corresponding element in each pannel. From the figure, it is 
clear that for ^^^Th no fragment is produced with isotopes beyond the drip-line. On the 
other hand, for ^^^Th a larger number of fragments are predicted to lie away from the 
drip-line, giving rise to some neutron emission simultaneously, which may be termed as 
multi-fragmentation fission. 
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